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We investigate the impurity resonant state induced by non-magnetic impurities in d-wave cuprate 
superconductors in two different impurity models: (i) in a pure potential scattering model within 
the T-matrix approach and (ii) in a Bose-Fermi Kondo model within the large- A/" formalism. We 
modify the superconducting host to resemble the nanoscale electronic inhomogeneities seen in the 
scanning tunneling microscopy (STM) spectra (i.e., the small- and large-gap domains) and study 
how this influences the impurity resonance. In the pure potential scattering model the resonant 
state is found to be a rather robust feature appearing in all host domains. On the other hand, in 
the Kondo-like model the impurity resonance can appear in the small-gap regions of the host and be 
completely suppressed in the large-gap regions; this is in agreement with the STM experimental data 
for Zn-doped Bi2Sr2CaCu2 0g+5. Our results imply that the Kondo spin dynamics of the impurity 
moment is the origin of the impurity resonant state in d-wave superconductors, rather than the pure 
potential scattering. 



I. INTRODUCTION 

Scanning tunneling microscopy (STM) experiments 
on surfaces of the high temperature cuprate supercon- 
ductor Bi2Sr2CaCu208+5 (BSCCO) provide a powerful 
tool for probing strong electronic correlations in this 
materialiii^'^'^ Since the STM measures a local density 
of states (LDOS) with a high real-space resolution it is 
widely used for studying impurity effects in c?-wave su- 
perconductors (SC).^ Extensive STM studies have been 
performed on both magnetic Ni and non-magnetic Zn 
impurities in BSCCO. Particularly interesting is the Zn 
resonant state, i.e., a large peak in the differential con- 
ductance (i.e., in the LDOS) at small but finite energy 
which is mainly localized near the Zn sitcjii^i^ Despite the 
fact that there are numerous theoretical studies dealing 
with the resonant state induced by non-magnetic impu- 
rities in d-wave cuprates, no consensus on the theoretical 
description has been reached yetn^^ 

The Zn^+ impurity has spin S = and it might be nat- 
ural to describe it as a delta-function potential scatterer 
I within a potential scattering (PS) modeljii^ In the PS 
model the resonance peak can be identified with a qua- 
sibound state whose energy can be tuned by varying the 
strength of the potential. However, a low-energy state 
appears only for large potential values, i.e., in the uni- 
tary limit. Moreover, some experimental findings cannot 
be easily explained within the PS model like the spatial 
dependence of the resonance peak which shows a global 
maximum on the Zn site and local maxima on the next- 
' nearest neighbor Cu sites 3 

Nuclear magnetic resonance (NMR) experiments are 
another powerful probe of strongly correlated systems 
like cuprates. In a series of NMR experiments it has 
been established that each non-magnetic impurity (e.g., 
Zn^+ or Li+) induces a net local 5=1/2 moment that 
is distributed on the Cu ions in its vicinityjiiii^ii^iii Fur- 
thermore, the effective induced moment interacts with 
the elementary excitations of the bulk material giving rise 
to Kondo-like physics which has indeed been observed in 



NMR experiments.^"* Consequently, one can treat a non- 
magnetic Zn impurity in the c?-wave SC within the Kondo 
model and identify the low-energy impurity state in the 
STM spectra with the Kondo resonance. ^^'^^ However, it 
is still not completely settled whether the potential scat- 
tering or the Kondo-like description is more appropriate 
for modeling the non-magnetic impurities in the cuprates 
and for explaining all experimental facts. ^ The main ob- 
jective of the present study is to resolve this question. 
In principle, Kondo-like behavior should be suppressed 
by an external magnetic field which would allow to iden- 
tify the origin of the observed impurity STM signal. So 
far, no such experiments have been performed and it has 
been shown that for this purpose rather high magnetic 
fields would be needed.J^ 

Recently the STM measurements on BSCCO have 
revealed the presence of nanoscale electronic inhomo- 
geneities which show up as spatial modulation of both the 
gap amplitude and the LDOS.^ii^ii^i^a These electronic 
inhomogeneities and variations of the LDOS attracted 
numerous theoretical and experimental studies 
Random out-of-plane dopant atoms which change the 
pair interaction locally^* and a competing order like 
antiferromagnetism'^'' have been proposed to be a pos- 
sible origin of the LDOS modulations in the cuprates. 
Importantly, the STM studies on BSCCO have identified 
two distinct electronic domains: (i) the domains with a 
small-gap magnitude and with well-developed SC coher- 
ence peaks and (ii) the large-gap domains without SC co- 
herence peaks. Remarkably, in the presence of impurities 
it was observed that impurity resonances are present only 
in the small-gap domains and no impurity resonances can 
be detected in the large-gap domainsi^ Using this exper- 
imental fact and motivated by the question about the 
origin of the impurity resonance in d-wave cuprate SC 
we investigate how the resonance state is influenced by 
putting the impurity into different domains of the SC 
host in both the PS model and in the Kondo model. 
To get quantitative results we employ a T-matrix and a 
large-A/" approach for the PS and the Kondo model, re- 
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spectively. We conclude that the Zn impurity resonance 
in BSCCO has predominantly a Kondo-like character. 

Very recently a new possibility for the origin of the Zn 
resonance peak in BSCCO has been proposed^i Namely, 
in d-wave SC short-ranged off-diagonal (t^) impurities 
generate low-energy Andreev states which are then as- 
sociated with the impurity resonance seen in the STM 
spectra. If the usual impurity scattering is absent, 
the Andreev resonant state exists symmetrically around 
the Fermi level. However, by including a weak poten- 
tial this symmetry is broken. This gives results which 
bear closer resemblance to the STM data. Below, we will 
comment on some of these results. 

The remainder of the paper is organized as follows. In 
Sec. HI] we introduce the relevant models for describing 
a non-magnetic impurity embedded in a d-wave super- 
conductor. In Sec. mil the T-matrix approach and the 
large-A/" formalism are explained. Our numerical results 
and a comparison with experimental data are presented 
in Sec. IIVI We conclude the paper in Sec. El 

II. THE MODEL 



B. The Impurity Hamiltonian 

We consider an effective impurity Hamiltonian which 
consists of a potential scattering and a magnetic term, 
^^imp = ^^pot + Hk- The potential scattering term de- 
scribes electrostatic interactions between the impurity 
atom and the conduction electrons whereas the magnetic 
term represents an impurity moment embedded in the 
SC host. 



1. Potential Scattering Model 

We assume that the potential scattering is com- 
pletely localized at the impurity site ro = (0,0), i.e., 
V{r) = VoS{r — rg). Consequently, the scattering oc- 
curs only in the isotropic s-wave channel and the scat- 
tering term is given by H^^^ = Vq J2a4)aCQa^ where 
Cq^ = N^^/^^^^Ckcr, and N is the number of sites. 
In the PS model the bulk Hamiltonian is defined as 
^^bulk = Hbgs- 



To describe a single impurity atom in a superconduct- 
ing material we will consider a Hamiltonian consisting of 
a bulk and an impurity part, H = -ffbuik + -ffimp- Both 
terms are specified below. 

A. The Superconducting Host 

As a model for the cuprate superconductor we use the 
usual BCS Hamiltonian 

^Bcs = E *k [(^k - f^)r' + \r^] (1) 

k 

where ^'j^ = (cj^^ "^-ki) ^ Nambu spinor with momen- 
tum k — {kx, ky), c\^^ (Ckcr) creates (annihilates) an elec- 
tron with momentum k and spin cr, and the three Pauli 
matrices in particle-hole space are denoted by t^'V^^. We 
use a tight-binding dispersion of the form 

Ej^ = — (cos kx + cos ky) + t2 cos kx cos ky + — (cos 2kx+ 

+ cos 2ky) + (cos 2kx cos ky + cos kx cos 2fcy) -I- 

+ cos 2/22; cos 2fcj,, (2) 

with values for the hopping elements that fit the 
ARPES data)^ h = -0.5951 eV, ^2 = 0.1636eV, h = 
-0.0519eV, = -0.1117eV, = O.OSlOeV, and fi = 
— 0.1305 eV for optimal doping {p = 19.5 %). The d-wave 
pairing function is given by Ak = Ao(cosfca; — cos ky)/ 2 
where Aq = 40meV. In order to model the nanoscale 
electronic inhomogeneities seen by the STM (i.e., the 
small- and large-gap domains) it is necessary to modify 
the pure BCS state; this will be discussed in Sec. IIV Al 



2. Bose-Fermi Kondo Model 

To specify the magnetic term in the impurity 
Hamiltonian we use the so-called Bose-Fermi Kondo 
modelj^i^liS^iS^ It describes a spin- 1/2 impurity, Sjmp, 
coupled both to the electrons of the conduction band and 
to a bosonic bath representing bosonic collective excita- 
tions in the host material. (Apart from its application to 
the impurity moments in the cupratesj^l the Bose-Fermi 
Kondo model has attracted considerable interest in the 
context of local quantum criticality within an extended 
dynamical mean-field theory,'^*' and also in certain meso- 
scopic systems like a quantum dot coupled to ferromag- 
netic leadsi^) 

In the d-wave cuprate SC it is reasonable to separate 
the low-energy degrees of freedom into fermionic Bogoli- 
ubov quasiparticles around the nodal points, and into 
bosonic antiferromagnetic spin-1 collective fluctuations 
near the ordering wave- vector Q = (tt, 7r)i^i^i^ The 
effective impurity moment induced by non-magnetic im- 
purity interacts with both the fermionic and the bosonic 
bulk excitations. Therefore, a non-magnetic impurity in 
cuprates is expected to be well described by the Bose- 
Fermi Kondo model^l 

^K = ^KSi,„p-So + 7oSin,p-0o- (3) 

Here, Sg = l/(2yV) Ekk'acr' cL'^<T<T'Ck'a' is the conduc- 
tion electron spin operator at the impurity site, <t are 
the Pauli matrices in spin space, and Jk is the usual 
Kondo coupling. The impurity moment is coupled to the 
bosonic bath via a coupling constant 70. The local field 
(pQ represents the local orientation of the antiferromag- 
netic order parameter and it is given by = J^qi^q-y + 
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b_^^)y^2JAq/ujq, whcre J is the bulk exchange constant, 
Aq is a dimensionless function containing the geometry of 
the system, and j — x,y, The simplest bosonic bath 
consists of free vector bosons, H^^^^ = J^q-y'^q^q-y^q-y^ 
with the dispersion Wq = + c^q^, where the momen- 
tum q is measured relative to the ordering wave-vector 
Q and c is the velocity of the spin excitations.^ At zero 
temperature, the effective bosonic mass m is equivalent 
to the bulk spin gap Ag (as seen in neutron scattering 
experiments^) . A spin gap Ag = corresponds to a bulk 
quantum-critical point controlling a transition between a 
quantum paramagnet (Ag > 0) and the antiferromagnet- 
ically ordered phase.— 



In the present study the Kondo Hamiltonian ^ de- 
scribes the physics of the so-called pseudogap (i.e., soft- 
gap) Bose-Fermi Kondo model. ^^'^^ There, a fermionic 
LDOS follows a power-law close to the Fermi energy, 
Pci^) ^ I^Ti where lj is measured from the Fermi level. 
Clearly, the rf-wave superconductor has r = 1. In addi- 
tion, at zero temperature and for Ag = 0, the bosonic 
LDOS also follows a power-law, ^ sgn(w)|w|^~^, 
where e = 3 — c?, and d is the dimension of the sys- 
tem (for cuprates d = 2). It is important to note 
that for 7o = the Bose-Fermi Kondo model ((21) re- 
duces to the extensively studied fermionic pseudogap 
Kondo model.— i^i^i^ Depending on the value of r and 
the presence or absence of particle-hole symmetry, the 
fermionic pseudogap Kondo model exhibits a quantum 
phase transition between a phase with Kondo screening 
for Jk > Jkc, and a phase where the magnetic moment 
is free (i.e., decoupled from the bath) for Jk < ^Kc- It 
is also known that the presence of an additional bosonic 
bath for 70 7^ 0, leads to a suppression of the Kondo 
screening."2^i21i^ Note that in the present analysis we will 
use the Bose-Fermi Kondo model with the bulk Hamil- 
tonian i/bulk = HbcS + Hbos^ 



To describe non-magnetic impurities like Zn (or Li) 
in the copper oxide planes it is appropriate to use the 
so-called extended magnetic impurity— NMR experi- 
ments have shown that the effective moment induced by 
a non-magnetic impurity is spatially distributed around 
the impurity site, Tq, and mainly resides on the four 
nearest neighbor copper sites.— Therefore, the Kondo 
model with an extended magnetic impurity is required 
in order to explain both the STM and NMR measure- 
ments on cuprates with non-magnetic impuritiesj ^^'^^ In 
the present calculation we consider a four-site Kondo 
model in d-wave superconductors where the effective 
magnetic moment is coupled to the four Cu sites s — 
(±1, 0), (0, ±1) adjacent to the impurity. One should 
bare in mind that in this four-site model the effec- 
tive LDOS of the fermionic bath seen by the impurity 
changes; this will be shortly discussed in Sec. IIII BT] 



III. METHOD 

We start with the definition of the conduction electron 
Green's function in the absence of impurities, which is 

given by £°(T,k) = -(r^*k(^)^k(0))- our nota- 
tion all underlined symbols denote matrices.) Its Fourier 

transform is given by g^°(w„,k) = J^'' dr (r, k) e*'^"'^, 
where f„ = (2n-|- 1)t:T represents the fermionic Matsub- 
ara frequencies and p is the inverse temperature. Note 
that the bare Green's function in real space is simply 
given by £°(iJ^„,r) = TV"! Eke*''"[«'^n " (ck - m)t" - 
Akr-]-i. 

A. T-matrix Approach 

The single scatterer problem, iJpot, can be solved ex- 
actly with the usual result for the perturbed Green's 
function^ 

£^(iz/„,r,r') =£°(ii^„,r-r') + 

+ gl{ii^n,r- ro)T(iz.„)£°(ii/„,ro - r'), (4) 
where the T-matrix is given by 

T{iiyn)^VoT'[l~Vogl{iiyn,ro)TY'- (5) 

The T-matrix depends on the strength of the scattering 
potential and on the local conduction electron Green's 
function, ^^(it'niro)- It is important to note that this 
local Green's function is the only input quantity repre- 
senting the bulk material; it differs significantly in the 
small-gap and the large-gap regions of the SC host (more 
details are given in Sec. IIV A"]) . 

B. Large- A/^ Formalism 

In order to obtain quantitative results in the Bose- 
Fermi Kondo model we use the large- A/" formalism. The 
large-A/" theory can be constructed by generalizing the 
spin symmetry of the impurity moment, the conduction 
electrons, and the vector bosons from SU(2) to SU(A/'). 
An antisymmetric representation of the impurity spin 
is given by S[^p = f^f^, - q^S^^, with the constraint 

fUf^ = 9o-^, where (Z^) creates (annihilates) aux- 
iliary fermions with /i = 1, . . . jA/" spin flavors. The con- 
straint is enforced by introducing a Lagrange multiplier 
Aq. The generalization of the spin symmetry leads to 
N"^ — 1 triplet bosons denoted by bq^f^'^ In what fol- 
lows we consider only the particle-hole symmetric case 
for which qq = 1/2. Therefore, in the large-A/" formalism 
the impurity Hamiltonian ^ is given by 

+ ^o(E4/m-9oA^)- (6) 
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It is important to emphasize that in the hmit A/" — !■ oo 
physics of both the fermionic Kondo model and the Bose 
Kondo model is controlled by a corresponding saddle 
point 

We proceed with defining the pseudo-fermion Green's 
function G'f'i'T) — — (T7-/(t)/^(0)) and its Fourier trans- 
form Q^f{ivn) = — Ao)~^- We drop the spin indices 
since the SU(Ar) symmetry is preserved. The full pseudo- 
fermion Green's function is determined from the Dyson 
equation 

gj\iun) = [e?(^^^„)]"' - S/(^^^„), (7) 

where the self-energy I]/(ji/„) has to be calculated. In 
addition, the Lagrange multiplier, Aq, can be determined 
by employing the occupation constraint for the pseudo- 
fermions, which yields 

i^go = (/V), where (/t/) ^ g^(r = 0"). (8) 

In order to determine E/(it'„) in the Bose- Fermi Kondo 
model we first briefly present the known results for two 
limiting cases, namely for (i) the fermionic Kondo, 70 = 
0, and (ii) the Bose Kondo model, Jk = 0. 

1. Fermionic Kondo Model: 70 = 

For 7o = by using a simple decoupling approximation 
one can obtain the non-interacting form of the Kondo 
Hamiltonian ®, i.e., = E^.[-("/mCo^ + ^.c.) + 

Ao/^/p] with the self-consistently determined condensa- 
tion amplitude v — Jx/J^'^^icQ^ffj}- This scheme, the 
so-called slave-boson mean-field approximation, becomes 
exact in the limit A/" — *■ 00 when the physics is domi- 
nated by a static saddle point. The known result for the 
pseudo-fermion self-energy is given bjs^^ 

EW(^l.„) =t;2^°(zz/„,ro), (9) 

where g^{iJ^n,ro) = Tr[£0 ro)(l + t^)/2] is the lo- 
cal Green's function of the fermionic bath; for Vq 
one should use ^ instead. The condensation amplitude, 
V, measures the interaction between the conduction elec- 
trons and the impurity spin. The non-zero value of v 
signals Kondo screening whereas for v = the moment 
is unscreened (free), v can be determined by calculat- 
ing the corresponding expectation value of the bosonic 
operator cj/ which is related to the mixed propagator 
gf,iT)^-(TrfiT)cl{0)). Hence, 

« = ^k(4/>> ^here (cj/) = ^/.(t = 0"). (10) 
A simple diagrammatic analysis leads to 

Gfdivn) ^ -vgf{iyn)Qc{i^n,'ro)- (11) 



The factor v arises due to the fact that no bare mixed 
propagator exists and the minus sign as a result of the 
anticommutation relations. A non-zero condensation am- 
plitude, w, at T = 0, decreases with increasing temper- 
ature and vanishes continuously at the so-called Kondo 
temperature, Tk. 

The equations (f7|)- (fTT|) are a self-consistent system of 
equations for determining the full /-propagator, the con- 
densation amplitude, w, and the Lagrange multiplier, 
Aq. These equations can be solved for a given Kondo 
coupling, Jk, and a local conduction Green's function 
at the impurity site. From the full /-propagator we 
can easily obtain the pseudo-fermion spectral function, 
Pf{Lo) = —lv[iQf{u) + i0+)/7r. A trivial solution of the 
self-consistent system of equations is given by w = and 
Aq = 0, which corresponds to an impurity spin which is 
completely decoupled from the conduction electrons. At 
high enough temperature this is the only solution. With 
decreasing temperature another non-trivial solution with 
V 7^ can appear. The non-trivial solution signaling the 
Kondo screening arises only below the Kondo tempera- 
ture and for Jk > Jkc- Note that Tk can be determined 
for a given Jk using PH)) in the hmit u ^ 0, Aq ^ 0; at 
zero temperature we use (jlO|) in the same limit to calcu- 
late the critical value of the Kondo coupling, Jkc- 

It is known that the slave-boson mean-field method 
has numerous artifacts at finite temperatures or in the 
vicinity of the quantum critical point for A/" = 2 which 
is the case of physical interest. Nevertheless, this ap- 
proach has been extensively used in the context of the 
pseudogap Kondo model^i^i^ and is known to capture 
the qualitative physics of the Kondo screened phase. 

We emphasize that in the four-site Kondo model, 
which we use in the present study, the LDOS seen 
by the impurity is given by the effective propagator 
Es,s' ¥'s¥'s'Tr[£^(^^/„, s, s')(1+t^)/2] . The Green's func- 
tion defined in Eq. ^ (i.e., if the potential scat- 
tering is present, Vq 7^ 0, we assume that it is lo- 
cated at the impurity site, ro) and Lps = +[— ]1 for 
s — ro = (±1, 0)[(0, ±1)]. For a more detailed discussion 
see Ref. il6|. 

2. Bose Kondo Model: Jk = 

For Jk = the Bose-Fermi Kondo model reduces to 
the so-called Bose Kondo model. In the large- A/" limit 
the pseudo-fermion self-energy is given by^ 

2 

where r»!^(r,ro) = -(T^(/)o(t)0J(O)) is the local Green's 

function of the bosonic bath. The self-energy Yi^^\ivn) 
is derived by summing all diagrams with non-crossing 
bosonic lines (the so-called "rainbow" diagrams) which 
is known as the non-crossing approximation (NCA) j^^i'^^ 
The self-consistency is obtained using (O, ([5]), and (I12p . 
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Interestingly, in the Bose Kondo model there is always 
a non-trivial interaction between the impurity and the 
bosonic bath, i.e., S]j'^(it'„) 7^ 0. At zero temperature 
and at the bulk quantum critical point {m — As — 
and ujq = cq) the self-consistent solution yields a power- 
law behavior for the pseudo-fermion spectral function, 
Pf{uj) ~ |a;|'^/^~^. This behavior does not correspond to 
the usual Kondo screening which occurs in the fermionic 
Kondo model, but rather reflects non-trivial fluctuations 
that are usually referred to as the bosonic fractional (fluc- 
tuating) phase. For finite but small Ag there are no bulk 
spin excitations at low energies. For u; ^ As the impu- 
rity moment is decoupled from the bath while the case 
w ^ As is analogous to the case Ag = 0; the crossover 
between these two regimes has been studied in detail in 
Ref. m. 



3. Bose-Fermi Kondo Model 

A generalization of the large- A/" approach for the Bose- 
Fermi Kondo model cannot be done in a straightforward 
manner and further approximations are necessary. (We 
note, however, that the large-A/" analysis has been ap- 
plied to a multi- channel version of the Bose-Fermi Kondo 
model42.) Guided by the fact that the physics of both the 
fermionic and the bosonic Kondo model is controlled by 
a saddle point in the A/" — > cxd limit, we combine both 
effects through the pseudo-fermion self-energy 

S/iz.„) = ^fiiiyn) + ^f\ivn)- (13) 

The fermionic part obtained within the slave-boson 
mean-field approach ([9]) and the bosonic NCA contri- 
bution (|f 2p are only formally decoupled. Their mu- 
tual effect can be easily understood through the self- 
consistency which we present below. Namely, diagrams 
containing both the conduction electron and the bosonic 
propagators are included through the full /-propagator in 
Y!"^\ivn)- In addition, according to Eq. (fTT|) the mixed 
propagator, Qfc, also contains the full pseudo-fermion 
propagator. Therefore, the interaction with the bosonic 
bath entering the mixed propagator through Qf influ- 
ences the condensation amplitude, i.e., Kondo screening. 

The system of equations d7|)-([T^ can be solved self- 
consistently in order to obtain the full /-propagator, w, 
and Ao for given couplings Jk and 70 , and for given local 
Green's functions of the fermionic and the bosonic bath, 
^°(w„,ro) and (ia;„, Tq), respectively; the results ob- 
tained by numerically solving this self-consistent system 
are presented Sec. IIVCI 

Before we show our numerical results in the next sec- 
tion, we discuss some limiting cases. For small coupling 
Jk and the fermionic bath exponent r > we expect no 
Kondo screening, w = 0, and that the problem reduces 
to the Bose Kondo model. On the other hand, for large 
Jk and small 70 the pseudo-fermion self-energy is dom- 
inated by its fermionic term and we should recover the 



results of the fermionic pseudogap Kondo model. Con- 
sequently, the two interactions of the impurity with the 
fermionic and the bosonic bath (i.e., the fermionic and 
the bosonic Kondo physics) compete; the same conclu- 
sion has been reached within the renormalization group 
approachi2Si22i2^ 

For m = As = there is a quantum phase transition 
between a fermionic Kondo-dominated phase, for Jk 3> 
7o, and a bosonic fluctuation-dominated phase, for Jk <C 
70- Note that there is no zero temperature free moment 
phase as long as 70 ^ and m = As = 0; this is the so- 
called bosonic fractional moment phase^^ An interesting 
quantity is Tk = ?!<;( Jk) where the Kondo screening sets 
in. It can be derived using Eq. (fTU)) provided that w — *■ 
and Ao — ^ 0. In our calculation we determine the critical 
value of the Kondo coupling, Jkc, in the Tk — > limit. 

IV. RESULTS 

In this section we present our numerical results which 
are obtained for a system at zero temperature. First 
we explain the way in which the LDOS of the host SC is 
modified. Afterwards, we show the results for both (i) the 
pure PS model and (ii) the Kondo model with and with- 
out the potential scattering at the impurity site. Finally, 
we compare them with the existing STM experimental 
data. 



A. Host Local Density of States 

As already noted recent STM measurements of the 
LDOS in cuprates have led to numerous studies; they 
deal with the possible origin of the nanoscale spatial 
variation"*'^^'^^ or identify the processes contributing to 
the tunneling spectra4i In the present calculation we do 
not attempt to explain the origin of the nanoscale inho- 
mogeneities but we model the host LDOS in a way to get 
close resemblance to the STM datai^i^ 

The STM experiments have established that the LDOS 
shows a nanoscale spatial modulation with two types of 
regions i^i^ (i) The first type corresponds to domains 
with small gaps and well-defined coherence peaks. In this 
region we assume that the BCS Hamiltonian ([T|) describes 
the bulk superconductor and we use the conduction elec- 
tron Green's function defined in Eq. (|4]) to obtain the 
LDOS. (ii) The second type of domains have large gaps 
and very broad gap-edge peaks whose amplitude is rela- 
tively low. In the following we will refer to these regions 
as the "pseudogaps" . 

In order to describe the host superconductor in the 
"pseudogap" domains we use a rather phenomenologi- 
cal approach. We write the conduction electron Green's 
function as [£^(w„,k)]"^ = [£°(w„,k)]"^ - E^(ii/„,k), 
where the bare propagator is the propagator in the pure 
BGS state. The main approximation enters through the 
form of the conduction electron self-energy, E^(w„,k), 
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which we do not calculate but we employ a simple model 
for it. (For a thorough discussion about the conduction 
electron self-energy we refer the reader to Ref. It 
is important to note that the self-energy is chosen in 
such a way that the corresponding LDOS obtained from 
the perturbed Green's function, Q_ (wnjk), mimics the 
LDOS measured by the STM on BSCCO,^i^i^ As al- 
ready mentioned, the local Green's function of the host 
SC is used further as an input (i.e., a given) quantity for 
the T-matrix and the large-A/" calculations. We assume 
that the imaginary part of the self-energy is given by 

S'/i(w,k) uJ^/Ao for \lu\ < Aq, and Sii(w,k) Aq 

otherwise!^ The w^-dependence ensures that the low- 
energy region of the LDOS (e.g., for uj <^ Aq) is not 
modified in all nanoscale domains which is in accor- 
dance with the STM measurements."^'^*' The real part 
of the self-energy is obtained by using the Kramers- 
Kronig relation. In addition, for simplicity we assume 
that Sii ~ S12 = S21 = ^22- 
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FIG. 1: (color online) The different local densities of states 
(LDOS) of the bulk superconductor which are used as in- 
puts for the T-matrix and the large-TV calculations. By going 
from a pure BCS state (denoted by "1") towards a "pseudo- 
gap" state (denoted by "6") the gap value increases and the 
superconducting coherence peaks disappear. The low-energy 
region remains unchanged. 

The host LDOS which we use in the present calcula- 
tion are shown in Fig. [1] and are denoted by "1, . . . , 6" . 
The case "1" corresponds to the pure BCS state with 
well-defined coherence peaks and a small gap while in 
other spectra the gap increases and the coherence peaks 
disappear. In this way we can model the STM spec- 
tra with small (large) gap regions with (without) SC co- 
herence peaks. Furthermore, the STM measurements 
show that the low-energy region close to the Fermi 
level remains unchangedj^ii^ this feature is also cap- 
tured in Fig. [TJ Note that larger gap values correspond 
to a lower hole doping level; we use the gap values 



(1.0, 1.1, 1.2, 1.4, 1.6, 1.8) Ao for the spectra "1, . . . , 6" , 
respectively. By adjusting the chemical potential /i one 
obtains the correct doping; more precisely we fixed fi = 
-(130.5, 125, 120, 115, 110, 105) meVfor the doping levels 
p =(19.5, 17.5, 15.9, 14.5, 13.1, 11.8) ^o.^'^ It is important 
to emphasize that the doping level does not significantly 
change our results and it can be in principle completely 
ignored. Nevertheless, we include the different dopings 
in order to model the experimental situation as close as 
possible. 



B. Potential Scattering 

We start our discussion with the pure PS model ( Jk — 
7o = 0). It is known that the delta-function potential 
scatterer generates a resonant state in the superconduc- 
tor LDOS and the resonance can be tuned towards the 
Fermi level by increasing Vb; the low-energy resonant 
state is possible only in unitary limit, i.e., when Vq is 
the largest energy scale in the problem. However, the 
PS model is not able to explain the spatial distribution 
of the resonance which is seen in STM experimentsi^i^ 
These results are summarized, for instance, in Fig. 1 of 
Ref. M- 

Here, we consider a non-magnetic impurity as a pure 
potential scatterer placed in different regions of the d- 
wave superconductor; those regions are characterized by 
the corresponding host LDOS (shown in Fig.[T|). In Fig. [2] 
we show the LDOS in a d-wave SC at the impurity site 
obtained within the T-matrix approach for a given Vq and 
different host LDOS. The potential scattering strength is 
chosen such that one gets a resonance in the low-energy 
region. 

Modifying the pure BCS state towards the "pseudo- 
gap" state does not change the resonance significantly. 
The resonance is slightly moved away from the Fermi 
level and its amplitude is somewhat suppressed. We con- 
clude that if the Zn resonance in BSCCO were of the PS 
type then it would be present in all regions (both small- 
and large-gap regions). This result is not surprising if we 
recall that in the T-matrix approach only the low-energy 
part of the host LDOS can influence the resonance close 
to the Fermi level. According to the STM experiments 
the low-energy part of the host LDOS remains unchanged 
(as in Fig. [T]) so that the Zn resonance in the PS model 
survives in all regions which is in disagreement with the 
STM data. Note that for larger Vq (i.e., for lower energy 
of the resonant state) the change in the position and the 
amplitude of the resonance state is even less pronounced. 



C. Kondo Response 

Here we present numerical results for the Bose-Fermi 
Kondo model obtained within the large-A/" approach. 
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FIG. 2: (color online) LDOS at the impurity site in the PS 
model for Vo = 0.2 eV and different host LDOS. A resonance 
appears at a low-energy scale ~ O(meV); it is slightly sup- 
pressed and moved away from the Fermi level by modifying 
the superconducting host, i.e., by going from "1" to "6" in 

Fig.m 



FIG. 3; A typical impurity spectrum in the Kondo screened 
phase for given fermionic and bosonic coupling constants and 
for a given spin gap. The scattering potential is Vo = 0.2, and 
the SC host is in the pure BCS state ( "1" ) . The impurity spec- 
tral function develops a peaked structure at a characteristic 
Kondo energy scale T*. 



1. Impurity Spectral Function 

In Fig. [3] the typical impurity spectrum in the Kondo 
screened phase {v ^ 0) is shown; it is obtained by solving 
the self-consistent system of equations ([71)-([T^ presented 
above. In the Kondo screened phase the impurity spec- 
tral function, pf{u!), develops a pronounced peak at small 
but finite energy, typically of the order of T* ~ C'(meV). 
Such a peaked feature does not develop in the local mo- 
ment phase where the Kondo screening is completely sup- 
pressed and V — 0. 

In the present large-A/" analysis Tk and T* are not 
equivalent, i.e., they do not behave in the same man- 
ner upon approaching the transition. Nevertheless, T* 
can be used to characterize Kondo screening since (i) the 
local magnetic susceptibility changes its behavior at T* 
and (ii) the peaked structure at T* also appears in the 
conduction electron T-matrix as seen both in the large-A/" 
and the NRG calculations J^ii^i^ (Note, however, that in 
the NRG calculations the characteristic energy scale, T*, 
can be identified with the Kondo temperature up to a nu- 
merical factor of order unityii^ii^) In the present work we 
perform a zero temperature calculation so that T* serves 
as a natural choice to characterize the Kondo screening. 
Furthermore, we associate the peaked structure in the 
impurity spectral function with the Zn resonance in the 
STM datai^ In other words, the Zn resonance is likely 
caused by a Kondo-like behavior of the induced effective 
magnetic moment with the Kondo temperature as a nat- 
ural low-energy scale. 

As already noted the resonance peak seen by STM 
does not decay monotonically with distance from the Zn 



site but rather oscillates, producing local minima and 
maxima in the DOS map. Remarkably, the four near- 
est neighbor Cu sites around the impurity have no DOS 
local maxima associated with them (in contrast to the 
next-nearest and next-next-nearest neighbor Cu sites) 
This non-monotonic spatial dependence cannot be natu- 
rally explained within the PS model. On the other hand, 
the four-site Kondo model for an impurity in a d-wave 
superconductor can capture this featurej ^^i^^ 

2. Characteristic Kondo Energy Scale 

Here we analyze how the characteristic energy scale 
changes with the Kondo coupling when the impurity is 
placed in different regions of the SC host. In Fig. 2] we 
show how T* depends on the Kondo coupling when a host 
LDOS is modified from the pure BCS to the "pseudogap" 
case (the corresponding LDOS are shown in Fig. [1]). We 
see that a larger value of the Kondo coupling is required 
in order to screen the impurity moment in the large-gap 
regions (i.e., the large-gap regions have larger critical 
Kondo coupling, Jkc, below which the moment is free 
and above which the moment is screened). In particular, 
for certain values of Jk the screening is possible only in 
the small-gap regions and not possible in the large-gap 
regions. In other words, the Kondo-like impurity reso- 
nance is suppressed by going from the pure BCS to the 
"pseudogap" state; this is an important result which is 
supported by the experimental data.— This fact cannot 
be explained within the PS model but naturally arises in 
the Kondo model. 

Note that the values of T* ~ 0{meV) also agree well 
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FIG. 4: (color online) The characteristic energy scale as a 
function of the Kondo coupling for different host regions (the 
corresponding LDOS are given in Fig.[T]). The potential scat- 
tering and the bosonic bath are absent: Vb = 0, 70 = 0. The 
pure BCS state (small-gap domain) has lowest Jkc- 



rity completely decoupled from both baths. Furthermore, 
fermionic and bosonic Kondo physics compete since the 
larger 70 values require larger Jkc in order to screen the 
moment. (Note that the same conclusion was reached 
within the RG approach for the Bose- Fermi Kondo model 
with As = .^*^1^^ ) Also, for larger spin-gap values the 
bosonic bath is less effective in destroying the fermionic 
Kondo screening. (This fact was proposed earlier to ex- 
plain the NMR data showing a suppression of the Kondo 
screening by underdoping.— ) 



200 



I 




^ 100 



with the Kondo temperatures estimated in NMR experi- 
ments ~ 0(10 K).^"' In Ref . [l6l it has already been shown 
that the presence of weak to moderate potential scatter- 
ing does not change the results of the four-site fermionic 
Kondo model. In contrast, large values of the scattering 
potential induce a large DOS on the neighboring sites of 
the Zn impurity. As a consequence the critical Kondo 
coupling is strongly reduced and the Kondo scale is in- 
creased up to temperatures ~ 0(100 K), which is in dis- 
agreement with the NMR dataJ"* On the other hand, 
the Kondo temperature is suppressed in the presence of 
the collective bosonic modes which are included in the 
Bose-Fermi Kondo model. In what follows we present 
the phase diagram of the Bose-Fermi Kondo model ob- 
tained using the large-A/" approach and also show the 
suppression of the Kondo screening in the presence of 
the collective bosonic modes. 



3. The Phase Diagram 

The projection of the phase diagram onto the Jk-7o 
plane is shown in Fig. [5] where the critical coupling Jkc 
is determined for different values of the spin gap and the 
host SC in the pure BCS state ("1"). For 70 = we 
recover a fermionic pseudogap Kondo model with a fi- 
nite value of Jkc which is, as expected, independent of 
the spin gap; below Jkc the impurity moment is free. 
For 7o ^ and Jk < Jkc we distinguish two different 
cases: (i) For the gapless bosonic bath, Ag = 0, the im- 
purity moment is not completely free but rather placed 
into a bosonic fractional (fluctuating) phasci^ (ii) A cou- 
pling to a gapped bosonic bath, Ag 7^ 0, leaves the impu- 
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FIG. 5: (color online) The projection of the phase diagram 
onto the Jk-7o plane; the critical coupling Jkc is determined 
for different spin-gap values and the host SC in the pure BCS 
state ("1"). For Ag = and 70 = (70 7^ 0) there is a tran- 
sition between the free (fractional) moment phase and the 
Kondo screened phase. For As 7^ there is only a transi- 
tion between the screened and the free moment phase. The 
potential scattering is absent, Vb = 0. 

In Fig. [6] we show the projection of the phase diagram 
onto the Jk-As plane. The critical Kondo coupling is 
determined for different values of the bosonic coupling. 
For As = (As 7^ 0) there is a transition between the 
fractional (free) moment phase for Jk < Jkc, and the 
Kondo screened phase for Jk > Jkc- By increasing the 
spin gap the number of low-energy bosonic excitations 
is reduced which makes the bosonic bath less effective 
in destroying the Kondo screening and as a consequence 
Jkc decreases. The bosonic spin gap is doping dependent 
and decreases with underdoping.'^— As a consequence, one 
expects that the bosonic bath stronger suppresses the 
Kondo screening in the region with lower hole doping 
which corresponds to the region with the larger-gap in 
LDOS. 

As mentioned above in the presence of the potential 
scattering the DOS on the neighboring sites of the Zn 
impurity increases which results in a decrease of the crit- 
ical Kondo coupling. For example, for 70 = 50 meV and 
As = 40 meV and the host SC in the region "1" , we have 
Jkc = (10L3, 100.4, 97.3) meV for Vq = (0, 0.2, 0.5) eV, 
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FIG. 6: (color online) The projection of the phase diagram 
onto the Jk-As plane; the critical coupling Jkc is determined 
for different values of the bosonic coupling and the host SC 
in the pure BCS state ("1"). For As = (As / 0) there is 
a transition between the fractional (free) moment phase and 
the Kondo screened phase. The potential scattering is absent: 
Vo = 0. 
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FIG. 7: (color online) The characteristic energy scale as a 
function of the average energy gap: Jk ~ 120 meV, Vb = 
0.5 eV, 70 ~ 25meV, As = 20meV (open squares); Jk = 
100 meV, Vo = 0.2 eV, 70 = 25meV, As = 40meV (solid 
circles) and As = 20meV (open circles); Jk ~ 100 meV, Vb = 
0, 70 = (solid squares); the lines are guides to the eye. By 
going from the small-gap region to the large-gap region the 
Kondo screening of the impurity moment vanishes. 



respectively. 



4- Suppression of the Kondo Resonance 

Here we discuss how the characteristic energy scale 
changes by modifying the superconducting host. Figure 
[7] shows the Kondo scale, T*, as a function of an aver- 
age energy gap Aq. The average energy gap is defined 
as Ao = (A_ -I- A+)/2 where the quantity A_ (A+) is 
the energy of the first peak below (above) the Fermi level 
in the host LDOS without impurity. For the pure BCS 
state Aq corresponds to a superconducting gap Aq while 
in the "pseudogap" state one can estimate it from the 
LDOS shown in Fig. [T] 

All sets of data points in Fig.[7]indicate that the Kondo 
effect is more suppressed in the regions with a larger en- 
ergy gap. Furthermore, the suppression of the Kondo 
screening is a monothonic function of the average energy 
gap and it is, in principle, possible to identify the maxi- 
mum value of the energy gap above which the impurity 
resonance cannot appear in the LDOS; this agrees with 
the STM measurements.'^ (Note that the amplitude of 
the Andreev resonant state in Ref. [l^ is a non-monotonic 
function of the energy gap; the largest resonance ampli- 
tude one obtains for intermediate gap values.) It is im- 
portant to note that the Kondo energy scale ~ O(meV) 
can be tuned to zero by simply modifying the host LDOS; 
we have shown that this is not possible in the PS model 
(Fig. [2). On the other hand, in the Kondo model the 
low-energy region of the bath LDOS is important for the 
Kondo effect, but through the self-consistent procedure 



also the LDOS away from the Fermi level influences the 
Kondo screening. Since the high-energy part of the bath 
LDOS differs significantly in small- and large-gap regions, 
it is natural to expect that the Kondo resonance will be 
affected when the host LDOS is modified. 

Figure [7] also suggests that the presence of the bosonic 
bath additionally suppresses the Kondo screening. For 
the spin gap, Ag, we use a constant value regardless of the 
specific host LDOS, i.e., regardless of doping level. We 
could, however, associate different spin gaps to different 
host LDOS since the spin gap is doping dependent!^ the 
underdoped regions (i.e., the regions with a larger energy 
gap Aq) have smaller Ag. In this case, the suppression 
of T* with Ao would be even more pronounced. 



V. CONCLUSIONS 

Motivated by the question about the origin of the im- 
purity resonance in d-wave SC we have compared two 
frequently discussed scenarios, i.e., the potential scatter- 
ing model and the Kondo model. We have investigated 
how the STM impurity resonant state is influenced when 
the SC host is chosen in such a way that it resembles the 
nanoscale electronic inhomogeneities seen by STM.— (i) 
If the impurity atom is assumed to be a pure potential 
scatterer we have shown within the T-matrix approach 
that the impurity resonance appears as a robust feature. 
The resonant state is rather insensitive to changes of 
the SC host. It is present in both small- and large-gap 
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regions which is in disagreement with the experimental 
data, (ii) The suppression of the impurity resonance in 
the large-gap domains observed by STM can be naturally 
explained by assuming that the impurity atom gives rise 
to Kondo-type physics. We have used the Bose-Fermi 
Kondo model and the large- A/" method to show that the 
characteristic Kondo energy scale is affected by changing 
the host properties (i.e., by changing the corresponding 
host LDOS). The impurity resonance resulting from the 
Kondo screening can appear in the host small-gap re- 
gions and can be completely suppressed in the large-gap 
regions; this is in agreement with the STM experimen- 
tal data for Zn-doped Bi2Sr2CaCu208+5. Therefore, we 
conclude that the Kondo spin dynamics of the impurity 



moment is the origin of the impurity resonant state and 
that the pure potential scattering model is not sufficient 
to capture the physics of the non-magnetic impurity in 
d-wave superconductors. 
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low-energy local DOS because the oj'^-dependence in the 
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